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In this paper we develop the new approach in time series analysis with a variable time step and present 
the results of quantitative and qualitative estimation of randomness and regularity, and the study of non- 
Markovian effects of the X-ray emission intensity of the microquasar GRS 19154-105. Our estimation is 
, based on the application of the theory of discrete non-Markovian stochastic processes and gives a wide 

set of quantitative characteristics and parameters of the studied system with a variable time step. This 
Q I set reflects the change of the effects of randomness and regularity, the alternation of Markovian and non- 

O ' Markovian processes in the initial time signals. Initially these characteristics were determined for the simple 

^ ^, model systems constructed by means of a molecular dynamics method. Further the calculation of these 

[/3 ' parameters was executed for the X-ray emission intensity of the microquasar GRS 1915+105. We have 

^ ', developed a theoretical scheme for the detailed analysis of the time series with a variable time step. The 

O ' new theory allows us to extract more detailed information about the X-ray emission of the microquasar GRS 

C/3 , 1915-1-105 from the experimental time series. The received set of the parameters allows to estimate accurately 

randomness and stochastic behavior, regularity and robustness, non-Markovian effects of statistical memory 
in the intensity of the X-ray emission. The analysis of the experimental data with a constant and variable 
time step is carried out by means of phase portraits, power spectra of the initial time and event correlation 
function (TCF and EOF), and memory functions (MF) of the junior order, the first three points of non- 
^ Markovity parameter. 

' Key words: discrete non-Markov processes, time series analysis, randomness and non-Markovity, X-ray 

^ ! intensity of GRS 1915+105 

(N| '. PACS numbers: 05.40.Ca; 89.75.-k; 95.75.Wx; 97.80.Jp 

O 
l> 

1 Introduction of the phenomenon of chaos in dynamic systems has 

Q ' allowed to take a new look at functioning of complex 

" ^ , The problems of randomness and stochastic behav- systems. Chaos is absence of order and it charac- 

^ ior, regularity and robustness have been in the focus terizes randomness and unpredictability of changes 

of attention in the studies of real complex systems in the behavior of a system and impossibility to de- 

of various nature over the past years. The analysis termine their origin and reasons. 
. !^ , of individual properties and characteristics of real While studying multiform manifestations of ran- 

^ ' complex systems is impossible without registration domness, the authors have defined chaos, and its the 

: and quantitative estimation of various information parameters describing chaotic and regular stochas- 

■ ■ ' on manifestation of chaosity and randomness. By tic processes in different ways. As a rule, the search 

the change of the measure of randomness or regular- for the most accurate criteria for the estimation of 

ity it is possible to judge about the complex dynam- randomness or regularities in the dynamics of real 
ics of the system and its evolution. The discovery 
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objects is carried out if we have the information on 
the behavior of the weh-known models of nonlinear 
dynamics such as, for example, Lorentz model, Pois- 
son model etc. Initially for the similar systems one 
calculates the parameters determining their chaotic 
dynamics. Further the given parameters are cal- 
culates for real physical objects. On the basis of 
the comparative analysis of these characteristics it 
is possible to come to certain conclusions about the 
evolution and further dynamics of the studied sys- 
tem. However the above mentioned models of non- 
linear dynamics do not carry sufficient information 
about the internal properties of real objects. The 
models similar to real physical objects allow to avoid 
these defects. The model description of a gas, liq- 
uid and solid by the methods of molecular dynamics 
simulations can serve as an example. 

In addition for the estimation of optimal qualita- 
tive and quantitative parameters and characteristics 
of non-Markovian effects and effects of randomness 
and regularity of X-ray emission intensity dynam- 
ics of the microquasar GRS 1915-1-105 we use the 
method of molecular dynamics. The galactic micro- 
quasar GRS 1915-1-105 was discovered as an X-ray 
transient in 1992 |T|, and has been observed to be 
extremely luminous ever since. This binary system 
contains a 14 Mq black hole |2] accreting from a 
late-type giant of mass 0.8 ± 0.5 Mq [3] via Roche 
lobe overflow. GRS 1915-1-105 is unique among ac- 
creting Galactic black holes spending much of its 
time at super-Eddington luminosities [4j. It is an 
extremely variable source, exhibiting dramatic, ape- 
riodic variability on a wide range of timescales, from 
milliseconds to months [5]. GRS 1915-1-105 is lo- 
cated on the Galactic plane at a distance of ~ 11—12 
kpc [21 [6] and suffers a large extinction of 25-30 mag 
in the visual band. The basic characteristics of the 
GRS 1915-1-105 binary system is the systemic veloc- 
ity which is 7 = —3 it 10 km/ the orbital period 
of the system is Torb = 33.5 ± 1.5 days [7]. Spec- 
troscopic observations in the near-infrared H and 
K bands identified absorption features from the at- 
mosphere of the companion (mass-donating star) in 
the GRS 1915-^105 binary. The detection of ^"^CO 
and ^^CO band heads plus a few metallic absorp- 
tion lines suggested a, K — M spectral type and lu- 



minosity class III (giant) [2]. Hard X-ray studies in 
20 — lOO/ceF band have shown erratic intensity vari- 
ations on time scales of days and months ^-[10]. In 
Ref. [11] the technique of differentiating and rescal- 
ing was applied to the GRS 1915-1-105 X-ray data. 
As a result the existence of a fundamental time-scale 
for the system in the range of 12-17 days was found. 
In Ref. [12] it was concluded, that microquasar 
GRS1915-I-105, as any other black hole system, may 
be chaotic in nature. One of proofs of the similar 
behavior of black hole systems is chaotic variabil- 
ity of the X-ray emission [13]. For the estimation 
of a degree of randomness or regularity of the X-ray 
emission of the microquasar GRS 1915-1-105 the cor- 
relation dimension of the system ^12j is often used. 

In this work we present a new method of re- 
search of randomness, regularity, robustness and 
non-Markovian effects in the X-ray emission inten- 
sity of the microquasar GRS 1915-1-105 on the ba- 
sis of the theory of discrete non-Markovian stochas- 
tic processes |14j-|16]. The effects of non-Markovity 
in real complex systems, natural |14t I15j . live |14j . 
[I6]-[I9], biological and physical [211 [23] 

are of special interest for the correlation analysis. 
This method allows to define the set of character- 
istics and parameters, which contain detailed infor- 
mation about non-Markovian effects and degrees of 
randomness of the X-ray emission of the researched 
system. Initially we have calculated the given char- 
acteristics and parameters for model systems by the 
method of molecular dynamics: for a low density 
gas with greater randomness and the stochastic be- 
havior in particle motion; for a high density gas; 
for a liquid near the triple point; for a solid that 
corresponds to a greater regularity in movement of 
particles. On the basis of the received results we 
have executed qualitative and quantitative estima- 
tion of non-Markovian effects and of stochastic and 
regular regimes in the initial signal with a constant 
and variable time step for the X-ray emission inten- 
sity of the microquasar GRS 1915-1-105. 



Nonlinear Phenomena in Complex Systems 



R.M. Yulmetyev et. al.: Non-Markov Statistical Effects of X-Ray Emission Intensity 



315 



2 The theoretical framework of 
the statistical theory of discrete 
non-Markovian stochastic pro- 
cesses 

In this section we present the basic concepts and 
definitions of the theory of discrete non-Markovian 
stochastic processes [Il]-[in], used here for the anal- 
ysis of a time series with a constant and variable 
time step. The theory is constructed on the discrete 
finite-difference presentation of the Zwanzig-Mori 
kinetic equations \24:\ [25] well-known in statistical 
physics. The theory of the discrete non-Markovian 
stochastic processes is widely applied to the analysis 
of real systems of physical, biological, live and social 
nature |14j-|19j.|23j. The dynamic, kinetic and re- 
laxation parameters and characteristics of this the- 
ory, contain detailed information about individual 
properties and qualities of the studied complex sys- 
tem. 



istration of the X-ray emission intensity, {N — l)r 
is the total time of registration of the signal, r is 
a discretization time step. In the researched time 
series the discretization time At = r = 1 day. The 
mean value of the variable {X), fluctuation 6xj, and 
dispersion o"^ can be presented as follows: 



7V-1 

{X) = j^^x{T + jt), (2) 
j=0 

5xj = Xj -{X), '^'^ = j^Yl 

j=0 



For the description of the dynamic properties of 
the studied complex system (the dynamics of corre- 
lations) it is convenient to use the normalized time 
correlation function (TCF): 



2.1 The basic concepts and definitions 
of the statistical theory of discrete 
non-Markovian stochastic processes 
for the analysis of time series with a 
constant time step 

As a rule, the time registration of any characteris- 
tics or parameters of a complex system is carried 
out at discrete time intervals of equal length. It al- 
lows to show the essential features of fluctuations 
and irregularities and also the features of various 
dynamic regimes in the initial time series. Some 
key definitions and concepts of the statistical the- 
ory of discrete non-Markovian stochastic processes 
for the analysis of the time series with a constant 
time step |14)-[16) are given below. 

Let us present the time dynamics of the X- 
ray emission intensity of the microquasar GRS 
1915+105 discrete time series xj of some vari- 
able X: 

X = {x{T),x{T + t),x{T + 2t),..., (1) 
x{T + kr), . . . ,x{T + {N - 1)t)}. 

Here T is the initial moment of the time the reg- 



^ N-l-m 

x6x{T + (j + m)r). 



where t = mr. Here 6xj,6xj^rn are fluctuations of 
the variable X at j,j + m step, correspondingly, 
is an absolute dispersion of the variable X. TCF in 
Eq. (3) satisfies the requirements of normalization 
and attenuation of correlations: 



lima(t) = 1, lim a{t) = 0. (4) 

t— »0 t— >oo 

It is necessary to note, that for a certain class of 
complex systems the requirement of correlations at- 
tenuation is not always possible. 

By the use of the Zwanzig-Mori projection oper- 
ators technique ^241 [25] it is possible to receive an 
interconnected chain of finite-difference equations of 
a non-Markovian type [ll]-[l6] for the initial TCF 
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\a{t) 



a{t) and the memory function Mi{t) {i = 1, 2, n): 

rra— 1 

: -rAi Mi{jT)a{t - jr) 

j=0 

+Xia{t), 

m—l 

-tA2 ^ M2{jT)Mi{t - jr) 

j=0 

+A2Mi(t), 

. . . , 



AMi(t) 
At 



j=0 



+AnM„„i(t), (5) 

wliere Aj are relaxation parameters, and parame- 
ters Aj form the eigenvalue spectrum of Liouville's 
quasioperator L: 

.< W„_iLW„_i > 



An ■ 
Ar, 



< |W„_i|2 > 

,< W„_iLW„ > 



(6) 



(7) 



< lW„_i|2 > 

Orthogonal dynamic variables W„ are received with 
the help of the Gram-Schmidt orthogonalization 
procedure: 

where 6n.m is Kronecker's symbol. To compare the 
relaxation time scales for the initial TCF a{t) and 
the ith order memory functions Mi{t) we use the 
non-Markovity statistical parameter e. Initially the 
given parameter was used for the analysis of the irre- 
versible phenomena in a condensed matter |26j-[29j. 
According to |14j-|16j we introduce the relaxation 
times of the initial TCF and the nth order memory 
functions: 

N-l 

Ta = At'Y a{tj 
j=0 
N-l 

TMi = At ^Ml(tj- 

j=0 



"J J, 



N-l 

TM„ = At ^ Mn{tj 
3=0 



(8) 



Then the spectrum of the non-Markovity parameter 
is defined set of dimensionless quantities 

{e*} = {£1,^2, En-l}, 
El = Ta/TM-^,£2 = TM1/TM21 • • • ) 



(9) 



Thus, the value e„ characterizes the comparison of 
relaxation times of the memory functions M„_i and 
M„. The non-Markovity parameter allows to divide 
all relaxation processes into Markov, quasi-Markov 
and non-Markov processes. The spectrum of the 
non-Markovity parameter defines the stochastic pe- 
culiarities of TCF. 

In work [14j the concept of the non-Markovity 
generalized parameter for frequency - dependent 
case was introduced: 



(10) 



where is a power spectrum of the ith order 

memory function: 



2 



Af-l 

/xo(i^) = |At a{tj)cos2TTut 

j=0 

. . . , 

N-l 

iJii{v) = I At A'Ii{tj)cos2'Kvtj^. 

j=0 

The above-mentioned equations (3)- (5) present 
the case of Zwanzig-Mori's statistical theory [24^ 
[53] for the discrete statistical complex systems. 
The statistical theory of discrete non-Markovian 
stochastic processes allows to reveal Markov and 
non-Markov effects, the effects of statistical mem- 
ory, effects of a dynamic alternation of stochastic 
and regular regimes in the initial time series for the 
X-ray emission intensity of the microquasar GRS 
1915+105. 

2.2 A new approach to the analysis of 
discrete time series with a variable 
time step 

In many real complex systems the registration of 
the initial time signal by different reasons is carried 
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out at time intervals of different length. To such 
systems we can refer the objects of an astrophysical 
and seismological nature, some biological and social 
systems, economic and ecological objects [301 ED 

In the given work we offer a new approach to 
the description of discrete non-Markovian stochas- 
tic processes with a time step of variable length. 
Such presentation of the initial time signal allows to 
find the dynamic development of the system, con- 
nected with a not real time scale but with its con- 
sistent presentation. The basic idea of this method 
consists in fixing individual events as a sequence of 
dynamic values. It allows to consider the dynam- 
ics of the system as a sequence of individual events. 
As an example, here we present the analysis of the 
time registration with a variable time step of the 
X-ray emission intensity of the microquasar GRS 
1915-^105. 

2.2.1 The basic concepts and definitions 
of the theory of discrete non-Markovian 
stochastic processes with a variable time step 

Let us consider the chaotic dynamics of the X-ray 
emission intensity as a sequence of events which are 
"non-uniform" on a time scale: 

^ = {6,e2,6,---,6,---,ejv}, (11) 

where the intervals of time Atij = ti — tj, j = i + 
1, j = i — 1, are unequal. Here presents the event 
at the moment ti which follows after the event 
i = 1, N is the number of the event. 

The mean value (E) , fluctuation , absolute dis- 
persion cj^ for the set of the random variable E are 
defined as follows: 



(E) 



1 ^ 



a 



1 ^ 

-y 



1 ^ 



(12) 



According to [30] -[32] we shall define the correlation 
dependence in a sequence of events (11) as follows: 

^ N-m 
"('^) = ~ X] ^^i^^i+rn- (13) 



The function introduced in the similar way a{t) 
presents the event (not time!) correlation function 
(ECF). The general requirements suggest, that EOF 
should have the properties of normalization and at- 
tenuation of correlations: 



lim a(n) 



1, lim a{n) 

n—foo 



0. 



(14) 



For the description of a discrete sequence of events 
we shall use Liouville's finite-difference equation of 
movement: 



Ae^(n) 
An 



(15) 



where ^j(n + 1) = U{n + l,n)^j(n), U{n + l,n) is 
an evolution operator. 

Thus, the left hand side of the Eq. (15) can be 
submitted as follows: 



^ = U{n+l,n)Ci{n)-Un) = 
An 

{U{n + l,n)-l}C,{n). 



(16) 



Let us present the set of values of the dynamic vari- 
able (5^i as a ^-component vector of the system's 
state: 

a) the vector of the initial state: 



(17a) 



b) the vector of the final state: 

{S^m, S^m+l,S£,m+2, SS,m+k}, (17b) 



where 1 < k < N . 

Using the standard expression for the scalar prod- 
uct of vectors and relations (13), (17a) and (17b), 
we receive the "event" correlation function (ECF) 
for a stationary sequence of events (fTT]) in the fol- 
lowing way: 



an 



(Ai(l)A^+„(n)) 



|Ai(l)|2) 



(18) 



2.2.2 Kinetic equations for discrete non- 
Markov processes with variable time step 

Let us write down the equation of motion (15) for 
the vector state: 



AA-^,(n) 
An 



iL(n, 1)A 



m+k 



n 



(19) 
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By the use of the projection operator technique we 
can spht Euchdean vector's space of state A(A;) into 
two mutually-orthogonal subspaces: 

A(A;) = A'{k) + A"{k), 
A'{k) = UA{k), A"{k) = PA{k). (20) 

Here projector operators 11 and P have the following 
properties: 

|A^(l))(Ai(l)I 



p = 1-n, n = n2, 

^ = P, UP = PU = 0. 



(21) 



It allows to split Liouville's equation (15) into two 
appropriate equations in two orthogonal subspaces: 



AA (n) 
An 

AA"(n) 
An 



iLiiA (n) + 2Li2A (n), (22a) 
: iL2iA'{n) + iL22A"{n). (22b) 



Here Lij = UiLUj presents the matrix elements of 
Liouville's quasioperator L = Ln + ^12 + ^21 + -^22- 
Solving the Eq. (22b) and substituting into the Eq. 
(22a), we shall receive: 



AA (n + m) 
An 



iLiiA (n + m) 



+iLi2{l + iAnL22rA (n) 

m 

-L12 J^{1 + iAnL22P AnL2iA'(n + [m - j]). (23) 
j=i 

Acting on the Eq. (23) by the operator of 
(A(1)|/(|A(1)|^) and taking into account an idem- 
potentity property of projection operators, we can 
receive the following finite-difference equation for 
the initial event correlation function: 



Aa(n) 
An 



iXia{n) - AnAi ^ Mi(j)o(n - j). (24) 



Assuming that An = 1, it is possible to formally 
solve this equation: 

m 

a(n + 1) = {iXi + l}a(n) - Ai Mi (j>(n - j) . 



Here Ai is an eigenvalue Liouville's quasiopera- 
tor. Relaxation parameter Ai and memory function 
Mi(j) are defined as follows: 



Al 



A 



(A^(l)^Ai.(l)) 
(Ai(l)Li2L2iAi(l)) 



Mi{j) 



(|AUl)P) 
(Ai(l)Li2{l + iAnL22VL2iAl{l)) 



{AUl)LuL2iAlil)) 

Equation (24) contains function Mi{j), for which it 
is possible to repeat the procedure submitted above 
and receive appropriate finite-difference equations 
of a non-Markovian type for memory functions of 
senior orders n > 1. To simplify the given proce- 
dure and generalize the received results one can use 
Gram-Schmidt orthogonalization procedure [33j : 

This operation allows to receive a new vector of 
state Ws, contained in the memory function Ms{j): 



Wo = Ai, Wi = {iL - Ai}Wo, 
W2 = {iL - A2}Wi - AiWo, . . . 



(25) 



For new orthogonal dynamic variables we re- 
ceive an interconnected chain of finite-difference 
equations of a non-Markovian type for the sth order 
normalized correlation functions: 



= -AnA2 V M2(j)Mi(n - j) 

/\r) ^ — ^ 

+iA2Mi(n), 



^^ = -AnA.+if;M.+i(j)M.( 



n-j) 



where 



Mi(n) 
Msin) 



+i\s+iMs{n), (26) 
(Wi(l)Wi(n)) 

(|Wl(l)|2) ' 

(W,(l)W,(n)) 



|W.(1)|2) 
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+1 



(|W.|2) ' 

• • • ; 

(W,[iLW,+i]) 
(|W.|2) • 



(27) 



The frequency - dependence of statistical spec- 
trum of the non-Markovity parameter for the case 
of time series with a variable time step will be de- 
fined as follows: 



/^i(i^) 



(28) 



where is a power spectrum for the ith order 

correlation function: 

f N Y 

Mi(^) = \ An Mi(n) cos(27rnz^) > , 

I n=l J 

. . . , 

(zy) = |An^Mi(n)cos(27rnz^)| . (29) 



3 The experimental data and the 
details of computer simulations 

We analyze here two types of the experimental data 
of the X-ray emission intensity of the microquasar 
GRS 1915+105 [M!- The first set of the data 
presents a one-day averaged time series of the X- 
ray emission intensity in the period from February, 
1, 1996 to September, 1, 2004 (a step discretization 
r = 1 day) . For the analysis of the given time series 
we use the statistical theory of discrete non-Markov 
processes with a constant time step (see Sect. 2.1). 

For the analysis of the non-equidistant time se- 
ries (the 2nd type of data) we use the statistical the- 
ory of discrete non-Markovian stochastic processes 
with a variable time step (see Sect. 2.2). As ex- 
perimental data [31] we use the time series of the 
X-ray emission intensity of the microquasar GRS 
1915+105 with a variable time step for the period 
from February, 1, 1996 to September, 1, 2004. To 
study the stochastic properties of the X-ray emis- 
sion we have carried out an additional study of four 



model systems using the method of molecular dy- 
namics. We have studied a model system (a g as, a 
liquid, and a solid) consisting of 2048 particles of ar- 
gon molecules. The particles interacted by Lennard- 
Jones potential V{r) = Ae{{a/r)^'^ — (cr/r)^} with 
parameters e/ks = 120K and a = 3. 405 A [35]. 
Here e is a well depth, cr is a interatomic distance. 
The simulation was carried out at constant temper- 
ature T* = 0.6 and at various densities n* = 0.1, 
n* = 0.5, n* = 1.0 and n* = 1.5. We made use 
of the "velocity Verlet" algorithm to integrate the 

14 



equations of motion [36] with a time step 10 sec. 



4 Randomness, regularity and 
non-Markov effects applied to 
the dynamic analysis of simple 
model systems and the X-ray 
emission intensity of the micro- 
quasar GRS 1915+105 

In this section we present a new method of quantita- 
tive estimation of randomness, regularity and non- 
Markov effects of a time series. Preliminary the 
calculation of the degree of manifestation of non- 
Markov effects and randomness or regularities in 
the dynamic movement of a particle in the given 
cell will be carried out on simple model systems: 
a low density g dense g liquid near triple 
point, a solid. The level of randomness and non- 
Markov effects is established for each model system 
and is carried out by means of a set of various char- 
acteristics. The set of data of characteristics and 
quantitative parameters contains reliable informa- 
tion about the degree of randomness or regularity 
and non-Markov effects in the researched model sys- 
tem. Then on the basis of the calculation of these 
characteristics we shall proceed to the analysis of a 
real process: the event variability of the X-ray emis- 
sion intensity of the microquasar GRS 1915+105. 
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4.1 Randomness and non-Markov effects 
in simple model systems 

The study of the dynamic features of behavior of 
real complex systems of a different nature in car- 
diology, neurophysiology, epidemiology, biophysics, 
seismology shows the existence of close connection 
between the first point of the non-Markovity param- 
eter and a quantitative measure of randomness or 
regularity of the measured signal. To establish this 
connection we present the results of the study of 
several simple physical models. The purpose of the 
similar study consists in detecting typical features of 
behavior of the non-Markovity parameter for model 
systems with a different degree of randomness or 
regularity. The given models were constructed by 
the method of molecular dynamics. As an exam- 
ple we have considered four Lennard-Jones model 
systems: a low density gas (T* = 0.6, n* = 0.1 in 
reduced units with parameters, where a = 3.405A 
is an interatomic distance, e/ks = 120K is the 
well depth; T = T*s/kB,n = n*/a^); a dense 
gas (T* = 0.6, n* = 0.5); a liquid near the triple 
point (T* = 0.6, n* = 1.0) and a solid (T* = 0.6, 
n* = 1.5). 

The time dependencies of the y-component veloc- 
ity for one particle in the studied cell is submitted 
in Fig. 1. On the basis of the comparative anal- 
ysis we can reveal a clear distinction between the 
behavior of particles in each model. For the case 
of a low density gas (see Fig. 1) a weak correlation 
between the particle velocity and time is observed. 
It is connected with great (in comparison with the 
size of particles) distances between particles and the 
collisionless regime of their behavior. Within a long- 
time limit weak correlations caused by interaction 
of two (three) particles are detected. In the given 
model the most chaotic behavior of a particle in all 
the studied systems can be observed. Thus, the y- 
component of a particle velocity of a dense gas (see 
Fig. lb) differs in strong correlations with time. 
The interval of fluctuation scattering is relatively 
fixed. The model shows "moderate chaotization" 
in the behavior of particles. The following model 
corresponds to a liquid near the triple point (see 
Fig. Ic) is characterized by the state of "moderate 
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FIG. 1. The time dependence of the y-component of 
particle velocity: a) for a low density gas (T* = 0.6, 
n* = 0.1, in reduced units); b) for a dense gas (T* = 0.6, 
n* = 0.5); c) for a liquid near the triple point (T* = 
0.6, n* = 1.0); d) for a solid (T* = 0.6, n* = 1.5). 
The weakest correlations correspond to the y-component 
of the particle velocity dynamics for a low density gas. 
This dynamics characterizes the most chaotic behavior 
of the particle motion in the studied models. With the 
increase of density of the model system the correlation 
in the behavior of particles become more intensive. A 
model of a solid corresponds to the greatest probability of 
interaction between particles and the strongest velocity 
fluctuations. 

regularity" in the motion of particles. In compar- 
ison with the previous models the correlation be- 
tween the motion of molecules amplifies noticeably. 
The amplification of correlations is connected to the 
increase of the density of the system, accordingly, 
to greater intensity of interaction of particles. The 
amplitude of fluctuations is found within a certain 
interval of values. The obvious regularity in the 
motion of real molecules corresponds to the model 
of a solid state (see Fig. Id). The diagram shows 
appreciable symmetry of fluctuations regarding the 
value characterizing the condition of equilibrium in 
the given model. A high degree of regularity is de- 
fined by a high degree of correlation of the states of 
particles. 

The phase portraits on plane projections of the 
two first orthogonal dynamic variables (Wq, Wi) are 
submitted in Fig. 2. Here various degrees of ran- 
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FIG. 2. The phase portraits on a plane projection of the 
two first orthogonal dynamic variables (Wq, Wi) for the 
considered model systems: a) for a low density gas; b) 
for a dense gas; c) for a liquid near the triple point; d) for 
a solid. The structure of phase clouds indicates a level of 
randomness or regularity of the studied objects. Higher 
symmetry and regular concentration of phase points con- 
cerning the center of coordinates are characteristic for 
objects with a high degree of regularity. Decreasing of 
concentration of the phase points near center of coordi- 
nates is connected with the increase of the level of ran- 
domness in the system behavior. 



domness and regularity in the movement of the par- 
ticles in the system is observed. In all the phase 
clouds various symmetry is revealed. In case of a 
low density gas the individual points of the phase 
portrait generate precise closed structures with a 
complex chaotic form. In the next models the phase 
points concentrate near to the center of the coordi- 
nate system. The density of the phase cloud grows 
with the growth of regularity of the model. Concen- 
trated orbits, located around the central nucleus, 
dissipate. The phase portraits of the solid model 
are characterized by the greatest concentration of 
the phase points near the origin of coordinates. It 
is necessary to note that with the transition of the 
system from a) to c) the fluctuation scale of variable 
W2 increases. The degree of changes of the orthog- 
onal dynamic variables noticeably differs owing to 
the difference between the degree of correlations of 
particles. 



FIG. 3. The power spectrum of the initial time corre- 
lation function /io('^) for four model systems. The pa- 
rameter a was calculated for each power spectrum. By 
estimating values a for each model system it is possible 
to reveal its connection with a level of randomness or 
regularity in the behavior of the system. 



In Fig. 3 the power spectrum of the initial time 
correlation function fiQ{i') for four model systems is 
submitted. The fractal parameter a of the power 
spectrum fio{i') ~ (where a < ) was calculated 
on all frequency scale for all model systems. For 
a low density gas this parameter amounts to a = 
— 1.82; for a gas a = —1.75. For a liquid it is equal 
to a = —1.63; for a solid a = —1.33. Thus, between 
the fractal parameter of the power spectrum no{i^) 
and the level of manifestation of a randomness and 
stochastic behavior (the parameter a grows with the 
increase of a regularity degree) there is appreciable 
interrelation. 

In Fig. 4 the frequency dependence of the first 
point of the non-Markovity parameter (further sim- 
ply the non-Markovity parameter), calculated by 
Eq. (10) is submitted. 

While analyzing the majority of natural [14t I15| 
and physical [23j systems earlier we found the fol- 
lowing feature in the behavior of the first point of 
the non-Markovity parameter. The value £i{0) ~ 
10^ -7- 10^ corresponds to the dynamic states of the 
systems with the greatest level of randomness and 
stochastic behavior. For the similar states pro- 
nounced Markov effects (manifestation of instanta- 
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FIG. 4. The frequency dependence of the first point of 
the non-Marlcovity parameter £i(j^) for the considered 
model systems. The analysis of various physical and nat- 
ural systems specifies a special value of the parameter 
ei(0) in the estimation of level of randomness or reg- 
ularity. The parameter £i(0) make it possible to judge 
about the level of randomness or regularity in the system. 
The maximal values of parameter ei(0) 10^ 10^) are 
characteristic for chaotic and fully randomized processes. 
Decreasing of the parameter £i(0) up to values ~ 10" re- 
fiects the increase of level of regularity. It is possible to 
come to the similar conclusions by estimating the values 
of this parameter for the submitted models. 



neous or short-range statistical memory effects) are 
characteristic. With the increase of regularity and 
robustness in the system the numerical value of the 
first point of the non-Markovity parameter on zero 
frequency decreases up to a unit ei(0) ~ lO". Non- 
Markov processes with amplifying effects of long- 
range statistical memory correspond to such states. 

The analysis of the model systems lead to the 
similar results. The parameter ei(0) for the model 
of a low density gas (see Fig. 4) constitutes 29.76. 
This is a maximal value of the non-Markovity pa- 
rameter characterizing a low density gas i.e. the 
model with the greatest level of randomness and 
strong Markovian effects. The value of this param- 
eter ei(0) for a dense gas equal to 17.77, for a liquid 
- 1.18, for a solid 1.14. With the increase of the den- 
sity of the system the numerical value of the non- 
Markovity parameter on zero frequency decreases to 
a unit (with the increase of the density of the system 



the non-Markovian effects increase). This testifies 
to the possibility to use the non-Markovity parame- 
ter as an informational measure of manifestation of 
randomness and regularity effects. 

Thus, the submitted model systems allow to de- 
fine the set of qualitative and quantitative char- 
acteristics for the analysis of a degree of random- 
ness and non-Markov effects in real complex sys- 
tems. These characteristics carry detailed infor- 
mation about randomness or regularity in the re- 
searched system. To these characteristics we can 
refer: the time correlation in the initial time series, 
the shape of the phase clouds {Wi = f{Wo)), the 
fractal parameter of the power spectrum fJ,o{i') and 
the first point of the non-Markovity parameter on 
zero frequency ei(0). The most authentic and in- 
formative parameter of the level of randomness and 
the manifestation of non-Markov effects is the sta- 
tistical non-Markovity parameter. 

4.2 The analysis of the X-ray emission 
intensity of the microquasar GRS 
1915+105 for the series with a con- 
stant time step 

The time series of the X-ray emission intensity of 
the microquasar GRS 1915-1-105 [M] during the pe- 
riod from February, 1, to September, 1, 2004 (the 
time discretization is equal to one day) is submitted 
in Fig. 5. Here we find the presence of quasiperiodic 
structures which are connected to a relative regular- 
ity of the signal within certain time intervals. It the 
end of the time series the quasiperiodic structures 
get most noisy therefore their form is destroyed. 

The phase clouds on the plane projections for 
various combinations of dynamical orthogonal vari- 
ables Wi, Wj (where i,j = 0..3) of the X-ray inten- 
sity of the microquasar GRS 1915-1-105 are shown 
in Fig. 6. It is possible to judge the character of the 
X-ray emission intensity by the shape of the phase 
clouds. The phase clouds are of an asymmetric type 
and consist of a centralized nucleus with a high con- 
centration of phase points and several points sur- 
rounding the nucleus. 

The power spectrum of the initial TCF /xo(i^) 
(Fig. 7a) and three memory functions of the ju- 
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FIG. 5. The time dependence of the X-ray emission in- 
tensity of the microquasar GRS 1915+105 with constant 
time step (step discretization t = 1 day). It is possible 
to discover quasiperiodic structures in the time depen- 
dence. The given structures bring relative regularity in 
the X-ray emission intensity. 



200 
100 


200 

<f~ 100 

-100 

-1( 

200 
tf" 100 

5 



W„[l] 



w ri]^"" 



200 r 
■f" 100 



200 r 

7~ 100 

-100 - 
-10( 

200 r 

100 
S" 



w ri]^"" 



100 
W, [x-'j 



FIG. 6. The phase portraits on the plane projections 
of the four first dynamic orthogonal variables of the X- 
ray emission intensity of the microquasar GRS 1915+105 
(for time series presented in Fig. 5). The phase clouds 
consist of the centralized nucleus with a high concentra- 
tion of phase points and individual points scattered on 
the perimeter. 



nior order //i(z^) (where i = 1,2,3) (Figs. 7b, c, d) 
for the intensity of the X-ray emission are depicted 
in Fig. 7. All the figures are submitted on a dou- 
ble logarithmic scale. The fractal parameter of the 



FIG. 7. The power spectrum of the initial TCF iioiv) 
(a) and three memory functions of the junior order (b, 
c, d) for the X-ray emission intensity of the microquasar 
GRS 1915+105 (for time series presented in Fig. 5). 
The fractal parameter of the power spectrum /Uo(z^) is 
equal to a = —1.41. This value is an intermediate value 
between the similar parameters for the models of a liquid 
(a = -1.63) and a solid (a = -1.33). 



power spectrum iJ,o{i') is equal to a = —1.41. This 
value corresponds to the intermediate quantity be- 
tween the similar parameters for the model of a liq- 
uid a = —1.63 and a solid a = —1.33. In the region 
of frequencies i/=5.7x 10~^ + 7.2x 10~'^ f.u.{lf.u. = 
1/r where r is a discretization time step) a series 
of dynamic peaks is detected in the power spectra. 
This frequency range corresponds to a time interval 
of r = 13.8 + 17.4 days. The maximal peak corre- 
sponds to frequency = 6 x 10"^ f.u., of r = 16.6 
days. 

The frequency dependence of the first three 
points of the non-Markovity parameter Sj, where 
i = 1,2,3 for the intensity of the X-ray emission 
of the microquasar GRS 1915+105 is presented in 
Fig. 8. Of special value is the first point of the 
non-Markovity parameter on zero frequency ei(0) = 
4.94. This value occupies an intermediate position 
between the appropriate values for a gas and a liq- 
uid. For all the frequency dependencies in the region 
of frequencies z/=0.057 + 0.072/.U. a series of peaks 
is detected. 

In Table 1 some kinetic (Ai, A2, A3) and relaxation 
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FIG. 8. The frequency dependence of the first three 
points of the non-Marlcovity parameter ei{iy) {i = 1, 2, 3) 
for the X-ray emission intensity of the microquasar GRS 
1915+105 (for time series presented in Fig. 5). The pa- 
rameter El (0) plays a special role in the quantification of 
chaotic and regular processes due to maximal informa- 
tion about the effects of Markovity and non-Markovity, 
short-range and long-range memory, randomness and 
regularity effects. For the X-ray emission intensity of 
the microquasar we have ei (0) = 4.94. 

parameters (Ai,A2,A3) for the X-ray intensity of 
GRS 1915-M05 with a constant time step are sub- 
mitted. Let us note, that the kinetic parameter |Ai| 
means relaxation rate of the studied system. Small 
values of kinetic and relaxation parameters signify 
the manifestation of stochastic effects, randomness 
and stochastic behavior in the registered time sig- 
nals of the X-ray emission intensity of the micro- 
quasar GRS 1915-M05. 



Table 1. Some kinetic and relaxation parameters 
(absolute values) for the X-ray intensity of GRS 
1915-1-105 (constant time step) 



1 Al [r-^] 


A2 


A3 [r-^] 




A2 


A3 [r-'] \ 


1 0.22 


1.20 


1.04 


0.10 


0.10 


0.03 1 



4.3 The analysis of the X-ray emission 
intensity of the microquasEir GRS 
1915+105 for a time series with a 
veiriable step 

The initial record of the X-ray emission intensity 
of the microquasar GRS 1915-1-105 as a sequence 
of events is submitted in Fig. 9. The given series 
essentially differs from the time series with a con- 
stant time step (the registration was carried out in 



the same period) by a visibly big set of the experi- 
mental data and consequently is more informative. 
It allows to define quantitative and qualitative pa- 
rameters and properties of the studied system with 
a higher degree of accuracy and reliability. 
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FIG. 9. The discrete sequence of events in the X-ray 
emission intensity of the microquasar GRS 1915-1-105 
(time series with variable time step, mean time dis- 
cretization is equal to 96 min). Each point in the given 
sequence represents a single event of emission. 

The power spectra of the initial EOF (Fig. 10) 
and three memory functions of junior orders (Figs. 
lOb-d) are shown in Fig. 10. All the figures are 
submitted on a double logarithmic scale. On the 
whole frequency region of the power spectrum /io(^) 
strong fractality with the exponent a = —1.53 (see 
Fig. 10a) is observed. Let us note that the frac- 
tal parameter for the series with a constant time 
step corresponds to a = —1.41. Hence the initial 
time signal registered with a variable time step is 
characterized by greater randomness and stochastic 
behavior in comparison with a signal with constant 
time step. Memory functions fii{iy), where i = 1, 2, 3 
(see Figs. lOb-d) manifest the similar fractal behav- 
ior in a wide frequency interval. 

In Fig. 11 the frequency dependencies of first 
three points of the non-Markovity parameter ei{iy) 
(where i = 1,2, 3) are submitted. The first point 
of the non-Markovity parameter Si{i') in the event 
presentation (see Fig. 11a) allows to reveal ad- 
ditional features in relaxation processes of the X- 
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FIG. 10. The power spectrum of the initial correlation 
function (a) and memory functions of the junior orders 
(b-d) of the X-ray emission intensity of the microquasar 
GRS 1915+105 for a non-equidistant time record. The 
power spectrum of the initial ECF for the event series is 
characterized by greater fractality than the power spec- 
trum of TCP for the case of a constant time step. In 
the power spectra of the memory function is observed 
the peak on frequency ly — 0.2/.U. It reflects also the 
comparative analysis of the fractal parameter for both 
cases. 



ray emission intensity of the microquasar GRS 
1915-1-105. The value of the parameter on zero fre- 
quency has increased more than twice and is equal 
to £i{0) = 11.65. The similar behavior of the non- 
Markovity parameter (for the series with a variable 
time step) is connected with the amplification of 
quasi-Markovian and stochastic effects in the X-ray 
emission processes. In Fig. 11 markovization effects 
are observed in frequency behavior ei(i^) peak 
on frequency u = 0.238 f.u (where If-u = 1/An, 
An =1 event). The values of the two subsequent 
points of the non-Markovity parameter (see Fig. 
lib, c) are close to a unit and are of identical kind 
for both time series. 

In Table 2 some kinetic (Ai,A2,A3) and relax- 
ation parameters (Ai,A2,A3) for the X-ray inten- 
sity of GRS 1915-1-105 with a variable time step are 
submitted. The parameter |Ai| it is equal to 0.15. 
Smaller rate of a relaxation for the time series with 
variable time step is connected to more obvious and 
more pronounced stochastic behavior. It may be, 



of 6 
4 
2 




e, (01-11,65 a) 


2.5 
2 




b) 


2 




c) 






2 






1.5 
























Oi 
W 1 






1 


J 








0.5 





0.5 





0.2 0.4 

v[1/An] 



0.2 0.4 

v[1/An] 



0.2 04 

v[1/An] 



FIG. 11. The frequency dependence of the first three 
points of the non-Markovity parameter of the X-ray emis- 
sion intensity of the microquasar GRS 1915+105 for an 
event correlation function (for series of events presented 
in Fig. 9). The first point of the non-Markovity param- 
eter ei(j^) (a) is characterized by the following features: 
the value on zero frequency is equal to ei(0) = 11.65; 
the peak on frequency v = 0.238 f.u. connected with the 
sharp amplification of Markovian (stochastic) effects is 
observed. 

that the system needs more time for return to a 
stable state of equilibrium. 

Table 2. Some kinetic and relaxation parameters 
(absolute values) for the X-ray intensity of GRS 
1915+105 (variable time step) 



1 Al [r-^] 


A2 [r-'] 


A3 


Al [r-'] 


A2 [r--] 


A3 [r--] 1 


1 0.15 


1.16 


1.03 


0.05 


0.12 


0.02 1 



Thus, the approach to the description of statisti- 
cal discrete processes in the event presentation sub- 
mitted in this work allows to received more detailed 
and clear picture of the stochastic processes occur- 
ring in complex nature systems. 

5 Discussion and Conclusion 

Energetic and spectral properties, quasi-periodic os- 
cillations (QPO) and the diverse temporal variabil- 
ity of GRS 1915+105 are the focus of numerous 
studies in last years. 

In the paper [37| the problem of the limits con- 
cerning the physical information that can be ex- 
tracted from the analysis of one or more time se- 
ries typical of astrophysical objects has been con- 
sidered. The work |38j has clearly established that 
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radio emission from GRS 1915+105 is intimately 
related to the presence of hard (power-law domi- 
nated) intervals in the X-ray light curves. This in 
turn physically implies a clear relation between a 
radiatively inefficient flow close to the black hole, 
and a synchrotron-emitting overflow or jet. This 
suggest that MHD effects could be responsible for 
the production and/or confinement of the jets found 
in system. For the first time in Ref. [39] were ob- 
served concurrently in GRS 1915+105 all of the fol- 
lowing properties: a strong steady optically thick 
radio emission corresponding to a powerful compact 
jet resolved with the VLBA, bright near-IR emis- 
sion, a strong QPO at 2.5Hz in the X-rays and a 
power law dominated spectrum without any cuttof 
in the 3 — AOOkeV range. 

J. Rodriguez et al. in the work |40j have 
presented the results of simulations INTEGRAL 
and RXTE observations of the microquasar GRS 
1915+105. They have focused on the analysis of the 
unique highly variable observation and have shown 
that they might have observed a new class of vari- 
ability. Then they have studied the energetic de- 
pendence of a low frequency QPO from steady ob- 
servations. 

A scenario for the variability of the microquasar 
GRS 1915+105 starts from previous works, lead- 
ing to the tentative identification of the accretion- 
ejection instability as the source of the low fre- 
quency QPO and other accreting sources [3T]. A 
model for the ~ 30 minute cycles often exhibited 
by GRS 1915+105 is determined by the advection 
of poloidal magnetic flux to the inner region of 
the disk, and its destruction by reconnection (lead- 
ing to relativistic ejections) with the magnetic flux 
trapped in the vicinity of the central source. This 
could be extrapolated further to understand the 
long-term variability of this and other microquasars. 

In the paper [42j authors discuss the possible ori- 
gin of the following behavior: the QPO spectra are 
well modelled with a cut-off power law except on 
one occasion where a single power law gives a satis- 
factory fit (with no cut-off at least up to ~ 40keV). 
The cut-off energy evolves significantly from one ob- 
servations to the other, from a value of ~ 21.8keV 
to ~ SOkeV in the other observations where it is 



detected. It was suggested in the work [32] that the 
compact jet detected in the radio contributes to the 
hard X-ray (> 20keV) mostly through synchrotron 
emission, whereas the X-ray emitted below 20keV 
would originate through inverse Compton scatter- 
ing. The dependence of the QPO amplitude on the 
energy can be understood if the modulation of the 
X-ray flux is contained in the Componized photons 
and not in the synchrotron ones. 

The variability pattern is characterized in Ref. 
[43j by a pulsing behavior, consisting of a main pulse 
and a shorter, softer, and smaller amplitude precur- 
sor pulse, on a timescale of 5 minutes in the JEM-X 
3 — 35keV light curve. It was revealed, that the 
rising phase is shorter and harder than the declin- 
ing phase, which is opposite to what has been ob- 
served in other otherwise similar variability classes 
in this source. The fit show the source to be in a soft 
state characterized by a strong disc component be- 
low ~ 6keV and Comptonization by both thermal 
and non-thermal electrons at higher energies. 

The source, which was observed 3 times in the 
plateau state, before and after a major radio and X- 
ray flare, showed strong steady optically thick radio 
emission corresponding to powerful compact jet re- 
solved in the radio emission corresponding to power- 
ful compact jet resolved in the radio images, bright 
near-infrared emission, a strong QPO at 2.5Hz in 
the X-rays and a powerful law dominated spectrum 
without cut-off in the 3 — 300A;ey range |44) . 

Relativistic jets are now in Ref. [45j believed to 
be a fairly ubiquitous property of accreting com- 
pact objects, and are intimately coupled with the 
accretion history. Associated with rapid changes in 
the accretion states of the binary systems, ejections 
of relativistic plasma can be observed at radio fre- 
quencies on timescale of weeks before becoming un- 
detectable. However, recent observations point to 
long-term effects of these ejecta on the interstellar 
medium with the formation of large scale relativistic 
jets around binary systems. 

In Ref. [36] authors by Rossi X-ray Timing Ex- 
plorer have found that as the radio emission be- 
comes brighter and optically thick, the frequency of 
a ubiquitous 0.5 — lOHz QPO decreases, the Fourier 
phase lags between hard (11.5 — 60keV) and soft 
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(2 — 4:.3keV) in the frequency range of 0.01 — lOHz 
cliange sign from negative to positive, the coher- 
ence between hard and soft photons at low frequen- 
cies decreases, and the relativistic amount of low- 
frequency power in hard photons compared to soft 
photons decreases. 

Energetic dependence of a low frequency QPO 
in GRS 1915+105 have been analyzed in the work 
[17] . The results presented could find an explana- 
tion in the context of the Accretion-Ejection Insta- 
bility, which could appear as a rotating spiral or 
hot point located in the disk, between its innermost 
edge and the co-rotation radius. 

In the present work the dynamic features in the 
behavior of a particle in the given cell in condensed 
matter and the X-ray emission intensity of the mi- 
croquasar GRS 1915+105 were considered jointly 
on the basis of the statistical theory of non-Markov 
processes. For this purpose we used two various ap- 
plications of the submitted theory as it allows to 
estimate chaotic and regular, random and stochas- 
tic, Markov and non-Markov processes defining the 
dynamics of the studied system. By taking into 
account the effects of discreteness, long-range and 
short-range memory, statistical relaxation now we 
are able to define a set of valuable parameters and 
characteristics, which contain detailed information 
about the properties of the studied systems con- 
nected with randomness. 

On the basis of the comparative analysis of the 
initial time series, phase portraits, power spectra 
of the initial TCF and the first point of the non- 
Markovity parameter we have found some typical 
features of the behavior of regular and chaotic, 
Markov and non-Markov components of stochas- 
tic processes in studied systems. The systems with 
high level of chaotization (a low density g dense 
gas) are characterized by weak correlation between 
the initial signals, the presence of complex chaotic 
closed structures in the phase portraits, the high 
value of the fractal parameter (according to its ab- 
solute value) of the power spectrum of TCF and the 
big value of the first point of the non-Markovity pa- 
rameter on zero frequency. In the similar systems 
the effects of short-range or instantaneous statisti- 
cal memory that correspond to Markovian processes 



appear. In regular stochastic processes strong cor- 
relations in the initial time signals, clear symmetric 
structures of phase clouds are observed. The last 
one consists of centralized nucleus of high concentra- 
tion. Here low values of the fractal parameter (ac- 
cording to its absolute value) of the power spectrum 
of the initial TCF, small values of the first point 
of the non-Markovity parameter on zero frequency 
are observed. In the similar systems the effects of 
long-range statistical memory can be seen, and re- 
laxation time scales of TCF and memory functions 
are comparable. The processes describing the given 
systems are non-Markov ones. 

The use of the above submitted technique allows 
to estimate non-Markov effects, robustness and reg- 
ularity, stochastic behavior and randomness of the 
X-ray emission intensity of the microquasar GRS 
1915+105 when the initial signal with constant time 
step is registered. Here the scale of time fluctu- 
ations, long-range effects of memory, discreteness 
of various processes and states, effects of dynamic 
alternation in the intensity of the X-ray emission 
play an important role. In the dynamics of the 
initial time signal of the X-ray emission intensity 
the fast change of various regimes, sharp and unex- 
pected alternation of various types of fluctuations 
and correlations were observed. Taking into account 
discreteness of the experimental data, statistical ef- 
fects of long-range memory and the constructive role 
of fluctuations and correlations we have received a 
detailed information about the properties and pa- 
rameters which characterize statistical properties of 
the fluctuating X-ray emission of microquasar GRS 
1915+105. 

Thus, in this work, a new approach to the descrip- 
tion of discrete non-Markov stochastic processes 
with a variable time step in the event representa- 
tion is presented. This approach is based on the 
consecutive use of the idea suggesting the existence 
of the event correlation functions. The similar cor- 
relation functions are new physical quantities de- 
termining probabilistic interrelation between a se- 
quence of events. When analyzing the time sig- 
nal (the X-ray emission intensity of the microquasar 
GRS 1915+105) with a variable time step the fol- 
lowing characteristics were determined: the power 
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spectra of the initial correlation function and mem- 
ory functions of the junior orders, first three points 
of the non-Markovity parameter, the fractal param- 
eter of the power spectrum of ECF and the value of 
the first point of the non-Markovity parameter on 
zero frequency. These dependencies and parameters 
allow to receive additional information about the 
properties and characteristics of the studied system: 
the quasi-Markov character of relaxation processes 
in the intensity of the X-ray emission, amplification 
of Markov effects on certain frequencies, the fractal 
character of the power spectra of ECF and memory 
functions. 

Finally, this paper only takes the first step in in- 
troducing the concept of event correlation analysis 
of time series and defining it in terms of quanti- 
ties that can be calculated from an experimental 
data. We believe, the method developed can form a 
basis to start formulating further meaningful ques- 
tions regarding the notions and presentations for 
real complex systems. 
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